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Abstract 

A field theoretic description for inclusive semileptonic B meson decays is formu- 
lated. We argue that large regions of the phase spaces for the decays are dominated 
by distances near the light cone. The light-cone dominance allows to incorporate 
nonperturbative QCD effects in a distribution function. A one-to-one correspon- 
dence with the heavy quark effective theory is developed, which can estimate the 
first two moments of the distribution function. These conditions are useful but not 
restrictive enough to specify the distribution function, which must still be determined 
from experiment. Several model-independent predictions, such as scaling, sum rules 
of the hadronic structure functions and relations among them, are made. General 
formulas for the differential decay rates on several variables are presented, which are 
used for calculating the electron energy spectra with an Ansatz for the light-cone 
distribution function. 



1 Introduction 



The study of semileptonic and inclusive B meson decays is interesting on several accounts. 
They are some of the simplest decays to study theoretically and experimentally. For 
semileptonic processes the structure of the lepton current is completely known and the in- 
clusive hadronic tensor involves a sum over all final states, which makes the decay products 
incoherent. Inclusive semileptonic B meson decays are useful for obtaining parameters of 
the electroweak theory. Two standard model parameters and V u \, are extracted from 
them. In fact, the charmless inclusive semileptonic decay of B mesons is the main ex- 
perimental source at present which determines V u b- This follows from a measurement of 
the endpoint spectrum of the charged lepton energy, whose better understanding is highly 
desirable. 

The hadronic tensor for the decays involves short and long distance contributions, where 
confinement effects are also important. There are, however, several features of the decay 
which make its theoretical analysis possible. The mass of the B-meson is large enough in 
comparison to 1 GeV, so that the mass of the virtual W boson is large enough to produce a 
commutator of two weak currents which are at light-cone distances relative to each other. 
This occurs for large but not all regions of phase space, which allows us to use the methods 
of deep inelastic scattering (DIS) and replace the commutator of the two currents with 
its singularity on the light cone times an operator bilocal in the quark fields. The matrix 
element of the bilocal operator between B-meson states contains nonperturbative QCD 
corrections. Its light-cone Fourier transformation is related to a distribution function, in 
direct analogy to DIS. The distribution function was discussed earlier jl], [3] and since it 
involves a heavy quark it is expected to be large in the region where the b quark carries a 
large fraction of the B-mesons momentum. 

The hadronic matrix element has another feature which was discussed recently. Namely, 
it contains a heavy quark and several properties can be established in the heavy quark 
effective theory (HQET) 0J]. In this paper we adopt the HQET and show that there is a 
correspondence between results obtained by other groups || |6], [7|, ||, |9|] based on the HQET 
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and moments of the distribution function defined on the light cone. The first two moments 
of the distribution function can be estimated, which determine the mean value and the 
variance of it. These two properties are useful but not restrictive enough to determine the 
distribution function completely. We find that in order to reproduce the decay spectrum, 
we need more detailed knowledge of the distribution function. For this reason we find 
it appropriate to introduce the distribution function over the whole range of its variable, 
instead of its first few moments, and try to determine it from experiment. 

In this paper we present a field theoretic prescription for inclusive semileptonic B meson 
decays and determine properties of the distribution function. We intend to justify, in terms 
of field theory, the parton model for inclusive B decays [[[], |2|, |3| and to identify system- 
atic procedures for improving upon the parton model predictions. The advantage of this 
approach is that model-independent aspects of the analysis can be clearly separated from 
model-dependent ones. Several model-independent predictions are obtained. Among them 
is scaling: away from the boundary of the phase space, the hadronic structure functions 
depend on the scaling variable 

£+ = [q ■ Pb + ^(q-P B ) 2 -M B (q 2 -m 2 q )]/M 2 B 

only, and not on the momentum transfer squared q 2 directly. 

The paper is planed as follows. In sect. 2 we give the kinematics and the general 
formalism for inclusive semileptonic B meson decays. In sect. 3 we argue that very large 
domains of phase space are dominated by light-cone distances of the two weak currents. 
The distribution function is introduced in sect. 4, while in sect. 5 the general properties 
and the physical implication of the distribution function are discussed. In sect. 6 we use the 
techniques of the operator product expansion (OPE) and the HQET to estimate moments 
of the light-cone distribution function. Sections 7 and 8 include predictions. Some model- 
independent predictions are presented in sect. 7. We give the triple differential decay rates 
for both B — > eu e X u and B — > eu e X c channels in sect. 8. In order to produce quantitative 
features of the data, we propose an Ansatz for the light-cone distribution function consistent 
with its known properties. We evaluate the electron energy spectra using the Ansatz. The 
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conclusions are in sect. 9. 



2 Kinematics 

We consider the inclusive semileptonic decays B — > WiX q , where I = e, fi, or r and X q is 
any possible hadronic final state containing a charm quark (q = c) or an up quark (q = u). 
The decays are produced by weak interactions. The partial decay width is given by 
i r/ 3 P r/ 3 P n r/ 3 P n 

^ - ^ w <£?k m tn 4^] - « - £*>■ w 

where g is the four-momentum of the virtual W boson, P and i£ denote the four-momentum 
and energy, respectively. Their subscripts B, I, is, and i correspond to the B meson, the 
charged lepton, the neutrino, and the final state particle, respectively. The summation on 
n in Eq.(l) and hereafter implies a sum over all possible final hadronic states \n). At the 
tree level, the decay amplitude is given by 

M = vJ^=u(Pl)-fO- - 7B)v{Pu)(n |j M (0)| B), (2) 

where V q b are the matrix elements of the CKM matrix and the relevant charged weak 
current is 

j^x) =: q(x)^(l - 7 5 )6(a;) : . (3) 

| B) is the B-meson state with momentum Pg and is normalized according to (B\B) = 
2Eb{2tt) 3 5 3 (0). Perturbative QCD corrections, which were studied in |H| [11], [1^| and more 



recently in |T3], [TJ]], are not included. They can be added in an analysis of the data as a 
perturbation. For unpolarized leptons the partial decay width can be written as 

G 2 F \V qb \ 2 d^cP-P, 
di ~ (2n)*E B L W ^2E l 2E U > (4) 

where L^ u is the leptonic tensor 

IT = 2(P/'P; + P£P» - g^P x ■ P v + ief" a pP?P!])- (5) 
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W^ v is the hadronic tensor 

WW = E / n tA^t (2tt) 3 5 4 (Pb - g - £Pi)<5 jt(0) n>(n |j M (0)| B>. (6) 

n J i=l l Z7r J ^-^i i=i 

It is useful to express the hadronic tensor in terms of a current commutator 

WW = / dV™<5 |^(y), jt(0)]| S>. (7) 

Thus the commutator of two weak currents is relevant in inclusive semileptonic B meson 
decays. Furthermore, we shall see that the commutator of currents [j ^(y) , near the 

light-cone y 2 = plays a central role in inclusive semileptonic B meson decays. 

Generally, the hadronic tensor can be decomposed by introducing the hadronic structure 
functions Wi(q 2 , q ■ Pb) with two scalar variables chosen to be q 2 and q • Pb, 

W, u (P B ,q) = -g,vW 1 (q 2 ,q-PB) + ^^W 2 (q 2 ,q-P B ) 



B 



ie^^W^q 2 , q-P B ) + ^W 4 (q 2 , q ■ P B ) 



M 2 °^ " 

+ PB " qU ^ PBU W 5 (q 2 ,q.P B ), (8) 

where M B is the B-meson mass. The hadronic structure functions Wi(q 2 , q-Ps) characterize 
the structure of the decaying B meson. Nonperturbative QCD effects for the inclusive 
process under consideration are incorporated in them. 

Finally, the triple differential decay rate is obtained from the kinematical analysis 

— = G F\ V i b \ 2 y»> w (g) 

dE l dq 2 dq Q 32^E B ^ {) 

where the contraction of the hadronic with the leptonic tensor yields 

ITWp, = 2(q 2 -M?)W 1 (q 2 ,q-P B ) 

+ [AP l ■ P B q ■ P B - 4(P Z • P B ) 2 - M 2 B q 2 + M 2 B M 2 ]W 2 (q 2 , q ■ P B )/M 2 B 
+2[(q 2 + M 2 )q ■ P B - 2q 2 P l ■ P B ]W 3 (q 2 , q ■ P B )/M 2 B 
+M 2 (q 2 -M 2 )W i (q 2 ,q-P B )/M 2 B 

+4M 2 (q -Pb-Pi- P B )W 5 (q 2 , q ■ P B )/M 2 B , (10) 



where M t denotes the charged lepton mass. There are three independent kinematical 
variables in this inclusive phenomenology, for which we choose Pi ■ Pb, q ■ Pb, and q 2 . 
We see that for the massless lepton case only three hadronic structure functions Wi(q 2 , q ■ 
PB),W 2 (q 2 ,q ■ Pb), and W 3 (q 2 ,q ■ P B ) contribute. We will proceed to investigate the 
hadronic structure functions in a way as suggested by the light-cone dominance. 



3 Light-Cone Dominance 

B mesons are heavy. This implies that its decay dynamics is analogous to that of deep 
inelastic scattering, that is, the light-cone dynamics dictates the inclusive semileptonic B 
meson decay. To see this we start with an analysis of the hadronic tensor. 

According to the causality requirement the commutator in Eq.(7) has to vanish for 
space-like y, i.e., \jn(y), jt(0)] = 0, for y 2 < 0, and hence the integrand in Eq.(7) has a 
support only for y 2 > 0. 

Taking q = (q , 0, 0, q^), we have 

q-y = qoyo - qm = + qa)(yo - ys) + ^(?o - <?s)(i/o + ys)- (n) 

The dominant contribution to the Fourier transform of the commutator in Eq.(7) comes 
from domains with less rapid oscillations, i.e. q ■ y — 0(1); hence 

1 . . 

yo - 2/3 : , (12) 

qo + q3 

l . . 

yo + y3 , (13) 

qo - q3 

and 

22222^22 1 1 ri A\ 

y = y - Vi - V2 - Vz < y - Vz ~ - 2 — 3 = ~o- ( 14 ) 

% ~ % q 

Therefore, the dominant contribution to the integral (7) results from the range < y 2 < 
l/q 2 . This implies that as long as q 2 is large enough, q 2 > q 2 , the decays take place near 
the light-cone y 2 = 0, where q 2 is a reference scale and experience shows that q 2 — 1 GeV 2 . 
The scale q 2 . should be determined ultimately by experiment. 



For inclusive semileptonic B meson decays q 2 varies in the physical range of 

Mf <q 2 < (M B - M Xmm )\ (15) 

where My is the minimum value of the final hadronic invariant mass. The light-cone 
domain q 2 . < q 2 < (M B — M Xmin ) 2 covers most of the phase space, since the B meson is so 
heavy that the interval (Mb — Mx mi „) 2 — q$ ^ c Lo~Mf. Subsequently, in inclusive semilep- 
tonic B meson decays the light-cone contribution dominates over all other nonperturbative 
QCD contributions. Contributions far from the light cone are suppressed dynamically and 
kinematically. The leading approximation of nonperturbative QCD effects should be more 
reliable for the charmless decays B — > luiX u , where Mx min is negligible, and/or for the 
decays to the final states containing a r lepton with the mass M T = 1.777 GeV . For 
both cases nonperturbative QCD contributions far from the light cone are more seriously 
suppressed kinematically. 



4 Distribution Functions 



Applying Wick's theorem one obtains 

\jp(y)Jt(x)] = qivh^i 1 - i5){b{y)X x )hu{^ - 75)5(3) 

-b{x) lv {\ - l5 ){q(x), q(y)h,(l - 75)6(2/)- (16) 

The quark pairs qq could be either charm or up quarks. In the case of charm quarks the 
matrix element between B mesons is very small, because there is no spectator charm quarks 
in B~(bu) and B°(bd) mesons. The up-quark bilocal operator has a matrix element only 
between B~ states, but this term is suppressed for the reason given after Eq.(23). For the 
matrix element we keep only the second term in Eq.(16). The hadronic matrix element of 
the current commutator becomes 

(B\[j,(y),jt(0)]\B) = -(B |5(0) 7 „(1 - 7 B ){?(0), q(y)h^ - 7 5 )%)| B). (17) 

After some calculation, the expression (17) is transformed into 

(B \j»(y)jl(0)} B) = 2{S, auP - ie mi/ p)[d a A q (y)]{B 1(0)^(1 - 7s )%) B), (18) 



where S mv p = g m g v p + g^g va - 9y, v g a p- \(y) is the Pauli- Jordan function for a free 
g-quark, 

i\(y) = J j^e-^efaW - m D- ( 19 ) 

In Eq.(18) the matrix element is separated in two factors. The first factor contains the 
light-cone contribution in the form of the propagator and the long-distance part is included 
in the reduced matrix element. The decomposition is Lorentz covariant and each factor 
can be calculated in the Lorentz frame of preference. This is analogous to deep inelastic 
scattering where the production of two currents at light-like distances is expanded in terms 
of operators times their Wilson coefficients, which are obtained from perturbative QCD. 

In the light-cone limit y 2 — > 0, the reduced matrix element in Eq.(18) can be expanded 
in powers of y 2 from the general arguments of Lorentz covariance and translation invariance: 

oo 

B}=4irP§Y.(y 2 TFn(yPB). (20) 

n=0 

We define next the Fourier transform of F n (y • P B ), 

MO =jd(y P B )e* y - PB F n (y ■ P B ). (21) 

Near the light cone only MO survives, defined as the Fourier transform of the reduced 
matrix element at light-like separations 

B) |j,2 =0 . (22) 

The physical implication and properties of this distribution function will be discussed in 
the next section. 

Finally, the leading contribution to the hadronic tensor is obtained from Eqs.(7), (18) 
and (22): 

W, u = 4(S, av/3 - ie^p) J dU{t)e{{,P BQ - q )S[(CP B ~ qf - m 2 q ](^P B - q) a P^, (23) 

where m q is the mass of the quark in the final state. Repeating the above steps for the 
first term in Eq.(16) we obtain a distribution function for the up quark, whose ^-variable is 
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(B 6(0)7^1-75)%) 



fit) = MO 



d(y ■ P B )e*y p ° (B 6(0)^(1 - 7s)6(y) 



either negative or very close to one. In the former case, the negative £ is outside the support 
of the distribution function and when £ is very close to one the distribution function for a 
light-quark is very small, so that this contribution will be neglected. 
Comparing Eq.(23) with Eq.(8), we are led to the expressions 

(£+,£-) = 2[/(£ + ) + /(£_)], (24) 
W 2 (Z + ,Z-) = :r-V[£ + /(£ + ) -£-/(£-)], (25) 

Ws{£ + ,Z-) = -/(£-)], (26) 

W 4 (£h£-) = 0, (27) 
^ 5 (£+,£-) = W 3 (£+,£-), (28) 

where the dimensionless variables £± are defined as 



£±=[q-P B ± j{q-P B ) 2 -Ml{q>-ml)]/Ml. (29) 

Hence, the light-cone dominance ascribes the hadronic structure functions to a single uni- 
versal light-cone distribution function. The variable £_ occurs for the first time in the 
decays of heavy particles and is a consequence of field theory. 

The triple differential decay rate (9) can then be written in terms of the light-cone 
distribution function 

+£ + [4P, • P B q ■ P B - 4(P, • P B ) 2 - M 2 B q 2 + M 2 B M 2 }/M 2 B 
-[(q 2 + 3M 2 )q ■ P B - 2(q 2 + M?)P t ■ P B ]/M 2 B 

-&.<->£-)!■ (30) 
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5 Properties of the Light-Cone Distribution Function 



We discuss here some important properties of the light-cone distribution function. The 
distribution function is normalized to unity: 

1 



/ = MP I mV ' P ^ V ' Pb ( B \KWb(1 ~ 7s)%)| B) |, 2=0 



2M| 



Pg(B 5(0) 7m (1-75)6(0) B) = l, 



(31) 



due to the conservation of the b quantum number. 

We consider next /(£) in the rest frame of the B meson. In this frame, 



/(0 = ^ / dy e* M ^(B |6t( )P L 6(y )| B), 



(32) 



where the left-handed projection operator P L — (1 — 7s)/2. We insert a complete set of 
states between quark fields, translate the yo dependence out of quark fields. Then we get 



/(0 = E <K M £ - £M fl - p°J | (m |6 £ (0) | B) 



(33) 



where 6_l = Ppo. So we see that f(£) obeys positivity. The state \m) is physical and must 
have < p° m < M B , thus f{£) = 0, for f < and f > 1. Therefore, the support of the 
light-cone distribution function reads < £ < 1. These results are valid in any frame due 
to Lorentz invariance, although they are deduced in the B rest frame. Furthermore, we 
observe from Eq.(33) that /(£) is the probability to find in the B meson a b quark with 
a momentum £Pb- This is the familiar probabilistic interpretation of the parton model, 
except it is written in the B rest frame rather than the infinite momentum frame. 

It will be convenient to expand the light-cone distribution function in terms of deriva- 
tives of delta functions, 

-l) n 



/(0 = E 



n=0 



-M n (0<^)(£-0. 



(34) 



Such an expansion is very singular and any finite number of terms cannot represent the 
differential decay width. The expansion is convenient for comparisons with operator prod- 
uct expansions which also generate sequences with singular terms. The coefficient M n (£) 
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is related to the nth moment about a point £ of the distribution function as follows 

It follows, now, that Mq(£) = 1, the mean value \x and the variance a 2 of the light-cone 
distribution function can be expressed by the moments: 

ti = M 1 (0)=i + M 1 (i), (36) 

a 2 = M 2 (ji) = M 2 {i) - Ml(i). (37) 

To sum up, our results so far are quite general. It is shown that the 6-quark distribution 
function inside the B meson introduced in the parton model []IL 0, BJ is related to the light- 
cone Fourier transformation of the bilocal operator between B-meson states (22). This 
makes clear the connection with the parton model in inclusive B decays. In the next 
section we employ the techniques of the operator product expansion and the heavy quark 
effective theory to estimate moments of the light-cone distribution function. 



6 Moments of the Light-Cone Distribution Function 

To estimate moments of the light-cone distribution function, we must calculate the hadronic 
matrix element (B 6(0)7^(1 — 75)6(2/) B), which involves long distances and hence brings 
in nonperturbative effects of QCD. The techniques of the OPE and the HQET provide a 
possibility to calculate it in a systematic way. 

We start with the light-cone OPE. Since the b quark is very heavy within the B meson 
we can extract the large mass scale 

b(y) = e- im ^b v {y), (38) 

where m;, is the b-quark mass and v is the velocity of the initial B meson, defined by 
Pg = Mbv^ and the rescaled b-quark field b v (y) is related to the effective field of the HQET 
by Eq.(43) below. Upon performing the light-cone OPE the hadronic matrix element 
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becomes 

(B 6(0) 7 ^(1 - 76 )6(y) 5) 



3 — imjt 



^E^r-^( B ^ya-Tb^^-'-^CO) 5) (39) 

n=0 n - 

with = iD^. S denotes the symmetrization. This OPE keeps the leading twist operators 
with higher twist effects being neglected. The advantage of this expansion is twofold. First, 
the Lorentz structure of the matrix element allows us to express it in terms of the B-meson 
momentum 

(B \b v (0)^(l - Ts)^^ 1 • • • k^]b v (0)\B) = 

2(C n0 P£P£ 1 •••P^ + E M 2 B C ni g^P£ ■ ■ ■ P^P^ ■ ■ ■ P£) + 

i=i 

terms with g^ . (40) 

Terms with can be omitted on the light cone. Second, we can estimate some terms 
in the HQET as shown below. Substituting Eqs.(39) and (40) into Eq.(22) we have 

/(0 = E MmE^)^- (4i) 

n=0 Hl i=0 M B 

Comparing with Eq.(34), the nth moment about the point £ = m b /M B of the light-cone 
distribution function is related to the expansion coefficients as following 

n 

M n (m b /M B )=J2C ni . (42) 

i=0 

It is straightforward to note that the moment M (rrib/MB) = Coo is exactly equal to 1. 

We now go to the second step further. We employ the HQET to estimate other expan- 
sion coefficients. In this effective theory the rescaled b-quark field b v (x) is expressed by the 
velocity-dependent heavy quark field h v (x) by means of an expansion in powers of l/m&, 

where is the covariant derivative. The effective Lagrangian is 



b v {x) = [1 + ^- + G(—)]h v (x), (43) 



jChqet = Kiv ■ Dh v + h v ^^-h v - h v gG ^ a — h v + 0(-^), (44) 
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where igG^ = [D^, D„]. The series expansion in powers of l/m b is now explicit. 

By virtue of the methods based on the HQET the expansion coefficients C n i in Eq.(40) 
can be expressed in terms of small quantities, proportional to powers to be 

precise, the order of C ni and hence that of the moment M n {m b / M B ) is expected to be 
(A-Qcr>/ m f>) n - Hence nonperturbative effects can, in principle, be calculated in a systematic 
manner. A few coefficients are calculated to be 



C21 



5m b 

3M B 
2m b 



E b + 0(A 3 QCD / 



ml 



B 



E b + 0(A 3 QCD /ml), 
—K b + 0(A 3 QCD /ml), 



3M, 
2m 2 b 



c 22 = 0, 



(45) 
(46) 

(47) 
(48) 



where the dimensionless parameters, which parametrize the nonperturbative phenomena, 
are defined as ITU 



K b 
G b 



2M, 



2M f 



in 



(B 



2m\ 



- gG aj3 a 



a/3 



4ml 



-h, 



B), 



(49) 
(50) 



with E b = K b + G b . Both parameters are expected to be order (AQ C D/m b ) 2 . 

According to Eqs. (45-48) and (42), the first two moments of the light-cone distribution 
function are 



M (m b /M B ) 
M^mt/Ms) 

M 2 (m b /M B ) 



1, 
m b 

w 



E b + 0(A 3 QCD /ml) 



B 



3M 



B 



-K b + 0{A% CD /mt) 



(51) 
(52) 

(53) 



Mo is exactly equal to 1. M\{m b jM B ) receives no contribution of order AQCD/m b . As 
a consequence, there are no nonperturbative QCD corrections to moments at the level 
A QC D/m b ; they arise first at order (A QCD /m b ) 2 . 
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Substituting Eqs.(52) and (53) into Eqs.(36) and (37), the mean \i and the variance a 2 
of the light-cone distribution function /(£) are estimated up to (Agc^/m;,) 2 corrections to 
be 

»='w B [1+Et) - (54) 

° 2 = (^n 2 f-E!)- (55) 

Therefore, the light-cone distribution function /(£) is sharply peaked around £ = /i ~ 
nib/Ms and its width is of order Aqcd/Mb, in agreement with rather general expectations. 

We can furthermore implement a numerical analysis. The parameter Gj, for the B 
meson can be related to observables M 

m b G b = ~(M B .-M B ). (56) 

The experiment determines Gb to be —0.0065. The parameter K b was estimated using a 
QCD sum rule flT5"| . It has a large uncertainty and its range could be Kb = 0.006 — 0.012. 
In Tables 1 and 2 we list /x and a 2 for various parameters mj and Kb, respectively. The 
numerical evaluation indicates that \x = 0.85 — 0.95 and a 2 = 0.003 — 0.007, should the 
b-quark mass and the parameter Kb vary within the limits 4.5 GeV < nib < 5.0 GeV and 
0.006 < K b < 0.012. 

A few remarks are in order: 

(1) The free-quark decay model is reproduced, if one keeps only the first term in the 
series (41), i.e., /(£) = 5(1;— nib/ 'Mb). The nonperturbative corrections manifest themselves 
in the second and higher terms in Eq.(41). 

(2) Important information about the light-cone distribution function is obtained by 
applying the OPE and the HQET. The first two moments of it are reduced to two accesible 
parameters Gb and Kb- However, it should be emphasized that the first few moments 
do not exhaust the information hidden in the distribution function, because they do not 
determine it completely. This point becomes obvious when one observes that any truncated 
resummation of the expansion in Eq.(34) cannot produce a smooth function. This is the 
origin of the singularity at the upper endpoint of the electron energy spectrum found in 
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H 0, [8| , where a truncated HQET-based OPE has been used. Moreover, a truncated series 
gives rise to a delta function in the hadronic tensor W^, which demands the decay to be 
described by quark kinematics instead of hadron kinematics. This brings about ambiguities 
particularly at the endpoint of the b — > u electron energy spectrum. In quark kinematics 
the endpoint of the b — > u electron energy spectrum lies at E e = mi,/2, while the actual 
endpoint should be E e = Mb/2 from kinematics at hadron level. 

As a matter of fact, an infinite number of terms in the light-cone OPE (39) must be 
included, and we cannot reduce our task to the calculation of a few matrix elements of 
lower dimension operators. It should be noted that, when more moments are taken into 
account, then higher dimensional operators of the OPE are involved; although qualitatively 
the moment M n {m b / M B ) is expected to be of order (AqcD/ m b) n in the framework of the 
OPE and the HQET, their hadronic matrix elements are much more difficult to calculate. 
We may conclude that the OPE and the HQET can serve as an useful technique for 
obtaining additional information on nonperturbative QCD, but they alone are not sufficient 
to determine the shape of the heavy quark light-cone distribution function. The limitations 
of this method must be complemented by other theoretical approaches. Alternatively, the 
required information may become available from experiment. For example, the light-cone 
distribution function can be determined from a measurement of the triple differential decay 
rate d 3 T/dE e dq 2 d(, + , as advocated in ||T6f . 

(3) A resummation of the operator product expansion has been performed recently in 
fl7| , |T8| , [19|| in order to eliminate the difficulties mentioned previously in the approach for 
inclusive B decays [||, [F], Efl. Their treatments using an operator product expansion in the 
context of the l/m^ expansion is different in an essential way from what we have formu- 
lated here. The distribution function is introduced in section 4 in a general way (without 
invoking the HQET), in analogy to DIS. In the present work the HQET is employed to 
estimate moments of the light-cone distribution function. To this end one may use other 
nonperturbative approaches (e.g. QCD sum rules). Moreover, the physical interpretation 
of distribution functions is distinct: the structure function f(k + ) defined in [19|] determines 
the probability to find a b quark with the light-cone residual momentum k + inside the B 
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meson, while the light-cone distribution function /(£) defined by Eq.(22) is the probability 
of finding a b quark carrying momentum £Pg within the B meson. The differences originate 
drastically different predictions. 



7 Model-independent Predictions 

The light-cone distribution function contains the long-distance physics associated with 
strong interactions of the b quark inside the B meson. Although we know some important 
properties of this function derived on general grounds and remarkable progress was made 
in calculating its first few moments, it cannot yet be determined completely from first 
principles. Nevertheless, it is interesting to draw model-independent results without relying 
on the quantitative aspect of distribution functions. 

It is convenient to define the scaling structure functions Fj 

F±(Z + ,q 2 ) = \w^ + , q 2 ), (57) 

F 2 (£+,g 2 ) = ^-±^W 2 ^ + ,q 2 ), (58) 

-W 3 (Z + ,q 2 ), (59) 

F^ + ,q 2 ) = ^-j^W^q 2 ). (60) 

The kinematical analysis shows that in the light-cone domain and away from the resonance 
region, namely in the kinematical region where our approach applies, /(£-) is expected 
to be relatively small and can be ignored, since the light-cone distribution function is 
sharply peaked around £ = fx « m b /M B , as established by the HQET analysis in the last 
section. We anticipated this, because the function /(£_) describes the creation of a quark- 
antiquark pair inside the B meson through Z-diagram and the virtual correction should be 
small. Then Eqs. (24-28) are simplified to 

Z + m-r,q 2 ) = m + ,q 2 ) = -Z + m + ,q 2 ) = -Z + F 5 (Z + ,q 2 ) = £+/(£+). (6i) 

Two important features of these expressions are: 

(i) the structure functions Fj satisfy scaling: they become functions of = [q • Pb + 
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4 







mq)}/Mg alone and are independent of the momentum transfer 



squared q 2 ; 



(ii) the structure functions are related to each other through the light-cone distribution 
function. 

Thus the structure functions Fj(£ + ,g 2 ) are measures of the momentum distribution of 
the b quark in the decaying B meson. The first result is the analogue of the Bjorken scaling 
in B-decays. The second one will be evidence for the spin-1/2 nature of charged partons 
(the quarks), i.e. the analogue of the Callan-Gross relation. Furthermore, using Eq.(61) 
and the normalization of the light-cone distribution function leads to the following sum 
rules: 



These results are very similar to those in DIS, since the behavior of the hadronic tensors 
for both cases is dictated by the light-cone dynamics. The remarkable thing is that these 
results follow in general without having information about the specific shape of the light- 
cone distribution function. It is important to keep in mind that the above results are valid 
up to perturbative and non-leading nonperturbative QCD corrections. As in DIS including 
the perturbative QCD corrections will lead to scaling violation: quantities which scale will 
be modified by powers of Inq 2 . In order to uncover these properties a detailed measurement 
of the differential decay rate d 3 T /dEidq 2 dq , Eq.(9), and hence of the structure functions 
Fi(^+,q 2 ) is essential. 



The results obtained so far take into account the lepton mass effects. Now we concentrate 
on the inclusive semileptonic B decay to electrons, in which the electron mass is negligible. 
The differential decay rates for B — > ev e X q in the B rest frame is then simplified from 



Jo Jo 4+ Jo 




(62) 



8 Electron Energy Spectra 
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Eq.(30) to be 

d 3 T _G 2 F \V qb \ 2 q -E e 



dE e dq 2 dq 4tt 3 ^/ q 2 + m 2 
where the variables £± given in (29) become 



{/« + )(2£«? + - q 2 /M B ) - « + ~ £_)}, (63) 



with q = c,u relevant to the b — > c and the 6 — > u decays, respectively. 

As already pointed out, appropriately including nonperturbative QCD effects allows us 
to use the actual kinematical limits at hadron level and to give the correct E e end point. 
They are given by: 

0< E e <^£(1-M|hl), (65) 

0< q 2 <2E e (M B - M ^-^ ), (66) 
q 2 q 2 + M 2 -M 2 

E - + k- * £ — fcr^- (67) 

In Figs. 1 and 2 the phase spaces for the b — > c and the 6 — > w decays are demonstrated, 
respectively, which show also domains of validity for our approach. In the resonance region 
bound-state effects in the final state become large. However, physical quantities, integrated 
over an appropriate phase space region, could be calculated reliably in our approach. 

In order to calculate decay distributions, one needs an Ansatz for the light-cone distri- 
bution function /(£) consistent with the general properties pointed out in section 5. We 
propose a parametrization for the light-cone distribution function with two parameters a 
and b as follows 

m^ J J^±- 2 , (68) 

where N is the normalization constant. In addition, constraints on the parameters a and b 
are imposed by the numerical evaluation of the mean value and the variance of the light- 
cone distribution function implemented in section 6 based on the techniques of the OPE 
and the HQET. Within the bounds of Tables 1 and 2, we find a = 0.002 — 0.016 and 
b = 0.86 — 0.97. Our light-cone distribution function is illustrated in Fig. 3. 
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Next, we use the distribution function (68) to compute the electron energy spectra for 
both b — > c and b — > u decays, shown in Figs. 4 and 5 respectively. We see that both spectra 
are smooth and go to zero at the endpoint. There is no sharp spike as E e approaches the 
upper endpoint. A desirable consequence is that perturbative QCD corrections on top 
of our spectra will give finite results, i.e. without endpoint singularities, because of the 
vanishing of the spectra at the endpoint. Therefore the endpoint behaviour of the spectra 
with perturbative QCD corrections are smooth for both b — > u as well as b — > c decays. 

Comparing with the free-quark decay model, we find that nonperturbative QCD cor- 
rections appear to be significant for the spectra in the endpoint region, referring to Figs. 4 
and 5. In particular, nonperturbative QCD effects become more important for the charm- 
less B — > eu e X u decay spectrum in the endpoint region. Our b — > u endpoint spectrum 
is considerably softer than the free-quark decay spectrum, which is finite and non-zero at 
the endpoint. 

9 Conclusions 

We studied the semileptonic and inclusive B meson decays using field theoretic methods, 
which justifies several steps of the parton model. We give in terms of a light-cone distribu- 
tion function the general formulas for the B — > lv\X q decays, keeping the mass of the final 
quark m q = m u or m c . These formulas should be useful for analysing the decay spectra 
and determining the CKM matrix elements \V c b\ and \V u b\. 

Additional properties of the distribution function are derived from the HQET. These 
are the first two moments. They are useful but do not determine the distribution function 
completely. The distribution function must be determined, presently, from detailed fits of 
the data. For this reason we discuss a new parametrization of the distribution function, 
which in the heavy quark limit reduces to a delta function and thus reproduces the free- 
quark decay model. 

There is now at our disposal a complete and reliable formalism for inclusive semilep- 
tonic B meson decays. It can be used directly to fit experimental data. This could be 
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done, preferably, by experiment groups, which can appropriately include their experimen- 
tal conditions. We are also working along the line trying to reproduce the electron energy 
spectrum and other features of the data. Early comparisons of the parton-model spectrum 
with experimental data are encouraging 0. 

The approach is more reliable for the charmless B — > ev e X u decays and shows that a 
large percentage of the events (more than 15%) has E e > 2.3 GeV, which may result only 
from the b — > u transition. It should be possible to analyse the endpoint spectrum in our 
approach including radiative corrections and obtain a reliable value for V u f,. 
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Table 1 The mean /x of the light-cone distribution function for several values of m b and 

K b . 



m b \ 


'a 


eV] 


4.5 


4.6 


4.7 


4.8 


4.9 


5.0 


K b 




0.006 


0.849 


0.867 


0.886 


0.905 


0.924 


0.943 


K b 




0.007 


0.849 


0.868 


0.887 


0.906 


0.925 


0.944 


K b 




0.008 


0.850 


0.869 


0.888 


0.907 


0.926 


0.945 


K b 




0.009 


0.851 


0.870 


0.889 


0.908 


0.927 


0.946 


K b 




0.010 


0.852 


0.871 


0.890 


0.909 


0.928 


0.947 


K b 




0.011 


0.853 


0.872 


0.891 


0.910 


0.929 


0.948 


K b 




0.012 


0.854 


0.873 


0.892 


0.911 


0.930 


0.948 



Table 2 The variance a 2 of the light-cone distribution function for several values of m b and 
K b . 



m b [GeV] 


4.5 


4.6 


4.7 


4.8 


4.9 


5.0 


K b 


= 0.006 


0.00288 


0.00301 


0.00314 


0.00328 


0.00342 


0.00356 


K b 


= 0.007 


0.00336 


0.00351 


0.00367 


0.00383 


0.00399 


0.00415 


K b 


= 0.008 


0.00384 


0.00401 


0.00419 


0.00437 


0.00456 


0.00474 


K b 


= 0.009 


0.00432 


0.00451 


0.00471 


0.00492 


0.00512 


0.00533 


K b 


= 0.010 


0.00480 


0.00501 


0.00523 


0.00546 


0.00569 


0.00592 


K b 


= 0.011 


0.00527 


0.00551 


0.00575 


0.00600 


0.00625 


0.00651 


K b 


= 0.012 


0.00574 


0.00600 


0.00627 


0.00654 


0.00681 


0.00709 
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Figure Captions 



1. The phase space for the b — > c inclusive semileptonic decay of B-mesons. The dashed 
and dotted lines encircles the domain of validity for our approach. The region between the 
solid and the dashed curves is the resonance region. 

2. Same as Fig.l, but for the b — > u decay. 

3. The light-cone distribution function (68). The parameters are taken to be a = 0.0076 
and b = 0.92. 

4. The electron energy spectrum in B — > ev e X c decays. The solid line is obtained in 
our approach for M Xmin = m c = 1.5 GeV, a = 0.0076 and b = 0.92. The dashed line 
corresponds to the free-quark decay model spectrum, using mj = 5.0 GeV and m c = 
1.7 GeV set by a fit to ARGUS data. 

5. Same as Fig.4, but for the b — > u decay. The solid line is predicted in our approach for 
M Xmin = m u = 0, a = 0.0076 and b = 0.92. The dashed line results from the free-quark 
decay model, using raj, = 5.0 GeV and m u = 0. 
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